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We study realistic models of relativistic radiating stars undergoing gravitational col- 
lapse which have vanishing Weyl tensor components. Previous investigations are gener- 
is- ■ alised by retaining the inherent nonlinearity at the boundary. We transform the boundary 
condition to an Abel equation of the first kind. A variety of nonlinear solutions are gen- 



erated all of which can be written explicitly. Several classes of infinite solutions exist. 
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1 Introduction 

The evolution of a radiating star undergoing gravitational collapse, in the context of general 
relativity, has occupied the attention of researchers in astrophysics in recent times. The deriva- 
tion of the junction conditions by Santos pQ has made it possible to obtain exact models of an 
interior spacetime with heat flux to match with the exterior Vaidya spacetime; at the boundary 
of the star the radial pressure is nonzero. A variety of exact solutions has been generated over 
the years to study the cosmic censorship hypothesis, gravitational collapse with dissipation, end 
state of superdense matter, dynamical stability of radiating matter and temperature profiles in 



*SSM and SDM dedicate this work to Peter Leach on his 65 th birthday 

^Permanent address: Department of Mathematics, Durban University of Technology, Steve Biko Campus, 
Durban, 4001 



1 



the context of irreversible thermodynamics. De Oliviera et al. [2] proposed a radiating model 
in which an initial static configuration leads to collapse. This approach may be adapted to 
study the end state of collapse as shown by Govender et al. [3]. Kolassis et al. [4j assumed 
that the fluid trajectories are geodesic and generated exact solutions. These assumptions lead 
to particular solutions which may be used to study the physical features of the model such as 
the relaxational effects of the collapsing fluid on the temperature profile in theories of causal 
thermodynamics [5]-[lU]. 

In a recent treatment Herrera et al. [H] proposed a model in which the form of the Weyl 
tensor was highlighted when studying radiative collapse. This approach has the advantage 
of simplifying the Einstein field equations. However, Herrera et al. were not able to solve 
the junction conditions; only an approximate solution was found. Maharaj and Govender [T2] 
showed that it is possible to solve the field equations and the junction conditions exactly. 
Their solution is expressible in terms of elementary functions and contains the Friedmann dust 
solution as a special case. It is interesting to note that Herrera et al. [13] showed that other 
classes of solutions in terms of the elementary functions are possible. The exact solutions 
in both P2] and [13] depend upon the introduction of a transformation that linearises the 
boundary condition. The purpose of this paper is to demonstrate that it is possible to obtain 
other models by transforming the boundary condition to an Abel's equation which is necessarily 
nonlinear. We explicitly find exact solutions to the Abel equation under particular assumptions 
and thereby demonstrate that conformally flat radiating stars contain a richer structure than 
previously suspected. 

The main objective of this paper is to show that we can generate radiating relativistic stellar 
models without having to eliminate the nonlinearity at the boundary. In Section 2, we describe 
the basic features of the model for a radiating star and present the relevant differential equations. 
Results generated in previous investigations are briefly discussed in Section 3. These have been 
obtained by introducing a transformation that leads to a linear equation at the boundary. In 
Section 4, we introduce a new transformation at the boundary that leads to an Abel's equation. 
We show explicitly that a variety of exact solutions can be generated from the Abel equation. 
Consequently a variety of new models for radiating relativistic stars, with vanishing Weyl 
stresses, are possible. The physical features of the solutions are briefly considered in Section 5. 



2 The Model 

We consider a spherically symmetric radiating star undergoing shear-free gravitational collapse. 
The line element for shear-free matter interior to the boundary of the radiating star is given by 

ds 2 = -A 2 dt 2 + B 2 [dr 2 + r 2 (d6 2 + sin 2 6d<j) 2 )} (1) 
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where A = A(t,r) and B = B(t,r) are the metric functions. The energy momentum tensor 
including radiation for the interior spacetime is 



Tab = (P + P)u a u b + pg ab + q a u b + q b u a 



(2) 



where the energy density p, the pressure p and the heat flow vector q a are measured relative to 
the timelike fluid 4-velocity u a = -jSq. The heat flow vector assumes the form q a = (0, q, 0, 0) 
since q a u a = for radially directed heat flow. 

The nonzero components of the Einstein field equations, for the line element ([T]) and the 
energy momentum ([2]), can be written as 
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The Weyl tensor has all components proportional to 



^2323 = -B 2 sin 2 9 
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which represent the tidal forces. For conformal flatness these components must all vanish so 
that C2323 = 0. This leads to a nonlinear partial differential equation which is easily solved so 
that 

A = (C^ty 2 + l)B. (4) 
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Now from (j3bl) and (l3c|) . and using we obtain 

B" B' 1 _ 
"B 7 ~~ ~B ~ r~ 

which is the condition of pressure isotropy. Equation (jSJ) is integrable and we get 
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where Ci(t),C2(t) and C%(t) are functions of time. The forms for the metric functions A and 
B given above generate an exact solution to the Einstein field equations (J3]). 

The interior spacetime (pQ) has to be matched across the boundary r = b to the exterior 
Vaidya spacetime 

2m(v 
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The hypersurface at the boundary is denoted by S. The junction conditions at S have the form 
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on matching p]) and 0. For our model the junction conditions jSJ) reduce to the 
nonlinear ordinary differential equation 
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(9) 



resulting from the (nonvanishing) pressure gradient across the hypersurface S. Equation (J9j) 
governs the evolution of a radiating star with vanishing Weyl stresses. To complete the descrip- 
tion of this radiating model we need to solve the remaining junction condition Q. 



3 Elementary solutions 

The governing equation (jHJ) is a highly nonlinear equation presenting a formidable mathematical 
task to solve exactly in general. Note that in previous attempts to integrate (Q assumptions 
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were made that effectively linearised this boundary condition. We briefly summarise the known 
results. 

Herrera et al. [H] assumed the following approximate forms for the temporal functions 



Ci = eci(t), C 2 = 0, C 3 



t 2 ' 
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where < e << 1 and a > 0, is a constant. With the assumptions contained in (jTOJ) , (J9J) yields 
the approximate solution 
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Note that on setting G\ = the solution reduces to a collapsing Friedmann dust sphere. 

Maharaj and Govender [12] were the first to determine a closed form solution for Q). They 
assumed that C% = C (a constant), C2 = and introduced the transformation C3 = u~ 2 so 
that ([9]) takes the linear form 

u - 2Cbu - (Cb 2 + l)Cu = 

in the new variable u. Three categories of closed form solutions in terms of the elementary 
functions were obtained depending on the nature of the roots of the characteristic equation. 

Herrera et al. [13] extended this treatment to obtain a wider class of solutions. They set 
Ci — C (a constant), Ci = aC^ and introduced the transformation C$(t) = u~ 2 {t) so that (Q 
can be written as 
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which is linear in u. Then equation (iTTj) admits three classes of solution given by 
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Case 3: (C - afb 2 + (Cb 2 + I) [a (a - 2C)b 2 + (C - 2a)} = 



C 3 {t) = (A + ^^exp 



-2{C-a)b 
ab 2 + 1 



where f3\ and 02 are constants of integration. These solutions reduce to the Maharaj and 
Govender [12] model when a = 0. Note that other transformations that linearise ([9]) are 
possible as indicated in [T3] . 



4 Abel Equation 

The nonlinearity and complexity in the boundary condition ([9]) is clearly evident. It is therefore 
remarkable that closed form solutions in terms of elementary functions have been shown to 
exist as shown in Section 3. These particular closed form solutions have been generated from 
linearised forms of the governing equation ([9]). A natural extension would be a study of the 
existence of nonlinear solutions to the differential equation (jHJ). Such classes of solutions, if 
they exist, are important in the study of nonlinear behaviour of the shear-free, conformally flat 
model. Consequently we seek classes of solutions which retain the inherent nonlinear structure 
of (Q. These have not been found in the past due to the inherent difficulties of coping with 
nonlinearity. 

Here we consider a particular nonlinear transformation which leads to exact solutions. It is 
convenient to replace the function C\(t) with 



U = C x b 2 + 1 

Then the governing equation (Q may be written with some rearrangement as 
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Equation ( Tl3l) is complicated, but has the generic structure 
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where we have set 
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The transformed equation (jl4jl is an Abel's equation of the first kind in the variable U. Abelian 
equations are difficult to solve in general. However, the advantage of utilising the transformation 
(j!2p is that (|14p is a first order differential equation in U. In the following we present a 
comprehensive mathematical treatment of ( fl4l) and derive several classes of solutions. 



4.1 Case 1:^ = 

The restriction ^4 = immediately gives 



C 2 b 2 + C 3 = a 



where a is a constant of integration. Then (fl3l) becomes 
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which is an algebraic equation in U. 

Two cases arise: U = or U 7^ in fll6p . We easily find: 
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This solution is particularly attractive since we have an infinite choice of C 3 and no integration 
is required. 



7 



4.2 Case 2: V = 

With T> = we have two possibilities: either 2C?b 2 — C3 = or C3 = 0. 



We firstly consider 2C 2 6 2 — C3 = 0. Then (|T3|) becomes 
l> + f/ 

This is a Bernoulli equation with solution 
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This is an infinite class of solutions depending on C3. 

Now we consider C3 = 0. The Abel equation ({TBI) becomes 
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where K' is a constant of integration. Therefore for the second case we have the solution 
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Again we have generated an infinite class of solutions depending on C 2 . 
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4.3 Case 3: C = 

Upon setting C = we obtain the equation 

G, 1 
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This equation is quadratic in C 2 which implies 



Co 



C 3 b±JC 3 b 2 -AC 3 (C 3 + C 3 b) 



2b 2 



Hence C 2 is a known quantity if the function C 3 is specified. 
The Abelian equation ffl3|) has the form 
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The equation is complicated, but may be written concisely as 
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The simpler equation ( 1201 has the form of a Bernoulli equation with solution 
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where K is a constant of integration. Consequently for this case we have the solution 

\ 



Ci 

c 2 
c 3 



1 

b 1 



e 2t /\C 2 b 2 + C 3 



(C 2 & 2 + C 3 ) 3 / 2 K" + ±f 



(C 2 b 2 +C 3 ) 4 



dt 



1/2 



/ 



C 3 b± yc 3 ft 2 -4C 3 (C 3 + C 3 fe) 
arbitrary function of time 



(22a) 

(22b) 
(22c) 



Again an infinite class of solutions is possible. 



4.4 Case 4: 

This is the most general case and corresponds to the situation for which all of the coefficients 
A, B, C and V are nonzero. Equation (|T4"|) can be written as 
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so that a variables separable equation is possible if ^, ^ and ^ are constants. Then the solution 
may be written as the quadrature 
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It is important to emphasize that the additional constraints generated by ^, ^ and ^ being 
constant simultaneously are not easy to simplify and in fact may not be consistent. Note that 
when B = and C2& 2 + C 3 = a we regain Case 1 with ^4 = and ( Tl3l) becomes a Bernoulli 
equation. Then solution (JTT|) is applicable. However, in general, when £> = 0, the quadrature 
fTSBl is applicable. 



5 Discussion 

Herrera et a/. [XT] obtained the equation ([9]) governing the gravitational behaviour of a radiating 
spherical star undergoing shear-free gravitational collapse by imposing conformal flatness to the 
model. Investigations of this model have thus far been confined to exact solutions of linearised 
forms of this equation. The nonlinear behaviour of relativistic stellar models is an inherent part 
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of realistic stars undergoing radiative gravitational collapse. A study of the physical features 
of these models hinges on the solution of the governing nonlinear equations. In this paper we 
have presented exact solutions of the governing equation in which the nonlinearity has been 
preserved. This has been effected by transforming the equation (Q into an Abel equation (JHJ). 
We have found several classes of exact solutions given in (ITT)) . (ITS)) . (IT9l) and (1221) retaining the 
nonlinearity of the model. Note that these generate an infinite family of solutions which allow 
for a systematic study of radiating relativistic spheres in different scenarios. 

It is important to observe that simple particular cases can be generated from our nonlinear 
models in Section 4. For example with „4 = and {/ / we may obtain the line element 

ds 2 = B 2 [-dt 2 + dr 2 + r 2 (d6 2 + sin 2 6d(f) 2 )] (24) 

where (1241) is in conformally flat form. Here for the simple case C3 = a in (ITT)) we obtain 
G\ = 0, C2 = and then B 2 is a constant; the Minkowski spacetime is regained. It is inter- 
esting to observe that the case C\ = in (TTTj) also arises when C3 takes the value 

_ 2a 2 f3e' 2t l b - a 
3 ~ 2a[3e- 2t / b -3 

where (3 is a constant of integration and C2 = (a — C^)/b 2 . Then we can write 

2 _ f b 2 {2a(3e- 2t ' b - 3) 1 2 
~ [ab 2 - 2ar 2 - 2b 2 a 2 (3e~ 2t / b _ 

This simple analytic form facilitates the analysis of the physical features of the model. We 
consider now some physical features which may be investigated in future work. With suitable 
choices of the arbitrary time functions, the luminosity radius 

L = (rB)z 

may be easily found. The quantity 

d hip 
dlnp 

gives a measure of the dynamical instability of the stellar configuration at any given instant 
in time. We can use this result to confirm that the centre of the star is more unstable than 
the outer regions. Of particular importance is the thermal evolution of the fluid. The causal 
transport equation in the absence of rotation and viscous stress is 

rh b q b + q a = -n(h b V b T + Tu a ) (25) 

where h ab = g a b + u a u b projects into the comoving rest space, Tis the local equilibrium temper- 
ature, k{> 0) is the thermal conductivity, and r(> 0) is the relaxation time-scale which gives 
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rise to the causal and stable behaviour of the theory. As shown in Maharaj and Govender [12] 
for a physically reasonable radiative stellar model (I25p becomes 

(3{qB)T-° + A(qB) = _ a T!^H (26) 

where A and B are the metric functions. Both the causal and noncausal solutions of (I26I) may 
be investigated in a simple model. 
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